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We study the asymptotic interaction between two half-quantized vortices in two-component Bose- 
Einstein condensates. When two vortices in different components are placed at distance 2R, the 
leading order of the force between them is found to be (log R/£ — 1/2) /R 3 , in contrast to 1/ J? between 
vortices placed in the same component. We derive it analytically using the Abrikosov ansatz and 
the profile functions of the vortices, confirmed numerically with the Gross-Pitaevskii model. We 
also find that the short-range cutoff of the inter-vortex potential linearly depends on the healing 
length. 
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I. INTRODUCTION 

Multicomponent condensations appear in many sys- 
tems in condensed matter physics and QCD, from 
multi-component or spinor Bose-Einstein condensates 
(BECs), superfluid 3 He, multi-gap superconductors to 
chiral phase transition or color superconductors in QCD 
at high temperature and/or high density. Especially, 
multicomponent and spinor BECs admit a rich variety of 
topological excitations: Domain walls 1], Abelian Q and 
non-Abelian Q vortices, monopoles jj, 2D Skyrmions 
[H, 3D Skyrmions @, % vortons [|[, knots @, and D- 
brane solitons [To| • See @, [TJ QjJ as a review. Among 
these topological excitations, quantized vortices in mul- 
ticomponent BECs are the most important subject, be- 
cause they are closely related to the problems not only 
in other condensed matter systems such as superconduc- 
tors, superfluids, magnetism, and liquid crystal, but also 
in electro- weak theory [13J , QCD and grand unified the- 
ories in high energy physics, neutron stars and cosmic 
strings in cosmology [3 Ell • 

Interactions between quantized vortices are important 
information to determine the equilibrium configuration 
and dynamics of many vortices. It is known that, in 
a single-component BEC, the asymptotic interaction en- 
ergy per unit length of two parallel vortex lines separated 
by a distance R is proportional to \og(L/R), where L is 
the size of the system [jq . Thus, the inter- vortex force 
has 1/R dependence. Vortices in a BEC resemble with 
global vortices in relativistic field theories A 
relation between them was studied in [20| where it was 
suggested that spinning global vortices on a Lorentz vio- 
lating background behave as superfluid vortices. Global 
vortices are regarded as global cosmic strings or axion 
strings in cosmology and the inter-vortex force between 
two global vortices was shown to be 1 jR [l8[ , coinciding 
with the one in vortices in a scalar BEC, scalar super- 
fluid, and the XY model. Global vortices also appear in 
QCD; in chiral phase transition of QCD at high temper- 



ature or high density |2l|, |22j or color superconductor of 
extremely high density QCD [23|, |24j|. Inter- vortex force 
at large distance R was derived analytically at the leading 
order as 1/R for color superconductor [2J|, and cosev/i? 
with a relative orientation a of two vortices in the in- 
ternal space for chiral phase transition [22| ; see [25| as a 
review. 

However, the analytic formula of the vortex- vortex 
interactions in multicomponent BECs are still missing. 
Two-component BECs are the simplest example of the 
multicomponent condensates and have also attracted 
much interest to study the novel phenomena not found in 
a single component BEC. Recent experiments provide a 
good ground of study on the vortex- vortex interaction in 
two-component BECs by tunings-wave scattering length 
via a Feshbach resonance [26-28J. The minimally quan- 
tized vortex in two-component BECs has the winding 
number one half of a singly-quantized vortex in scalar 
BECs, and thus is often called a half-quantized vortex. 
Its mass circulation is fractionally quantized when mass 
densities of two condensates are different. Such a quan- 
tized vortex in two-component BECs has a composite 
structure, where a vortex core in one component is filled 
by the density of the other component. This vortex struc- 
ture was created experimentally through coherent inter- 
conversion between two components [29]. Interactions 
between the vortices in the different components are non- 
trivial because the two components interact only through 
the density, so that the vortex winding around one com- 
ponent does not directly experience the circulation of the 
other vortex winding around the other component. This 
fact results in an indirect interaction, where the filling 
component of each vortex core is affected by the circula- 
tion created by the vortex in the same component, drag- 
ging the vortex in which it is filled. Although the inter- 
active dynamics of two vortices in two-component BECs 
was studied numerically by Ohberg and Santos [3(j, the 
analytical form of the interaction force was not discussed. 

In this paper, we consider the asymptotic interaction 
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between two vortices in two-component BECs. We con- 
sider the vortex-vortex interaction for two cases (i) two 
vortices are placed in the different components and (ii) 
those in the same component. For the case (i), the lead- 
ing order of the inter- vortex force between them at dis- 
tance 2R is found to be (\ogR/£_ - l/2)/R 3 with the 
short-range cutoff £, in contrast to the one 1/R for the 
case (ii) and vortices in a single-component BEC. We de- 
rive it analytically using the Abrikosov ansatz and the 
asymptotic profile functions of the vortices. We then 
confirm it numerically. We also find that the short-range 
cutoff £ of the inter- vortex potential linearly depends on 
the healing length. 

This paper is organized as follows. Sec. II is devoted 
for deriving the analytic form of the asymptotic inter- 
vortex force. In Sec. Ill we confirm the analytic results 
obtained in Sec. II by numerical calculations of the Gross- 
Pitaevskii equation. Summary and discussions are in 
Sec. IV. In Appendix, we describe some details of the 
calculation of integrals in Sec. II C. 



II. STATIC INTER- VORTEX FORCES 



A. The model 



We start with an energy functional for two-component 
BEC system 
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where ^ is a condensate wave function of the i-th com- 
ponent (i = 1,2) with mass m,. The coupling constants 
gi,g 2 and g\ 2 stand for the atom-atom interactions; the 
*i and ^2 components repel or attract for g i2 > 
or gi2 < 0, respectively. The coupled Gross-Pitaevskii 
(GP) equations are obtained by the variational principle 
ihd t ^ l = 5E/5^* as 
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The stationary coupled GP equation is given by consid- 
ering a time dependence ^(x, t) = e _v<i/,ft *j(x) with 
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The potential energy V with the quadratic terms 
— /xi|$i| 2 — /J-2 1 ^2 1 2 induced by the chemical potential is a 
quadratic function of X = | 2 > and Y = \^ 2 \ 2 > 0, 



V(X,Y) = ^X 2 
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g 12 XY - Ml X - M2 Y. (8) 



Let gi,g 2 be positive, then the potential V has a mini- 
mum when 



A = VxxVyy - Vxy = 9\9i - 3i2 > °; 
^\9i ~ ^i9\2 > 0, /u 2 3i - ^g 12 > 0. 
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The amplitudes of the ground state arc then given by 
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In the following, we consider the situation in which the 
above inequalities are satisfied. Since there are two con- 
densates, two U (1) symmetries are spontaneously broken. 
Accordingly the order parameter space is 



T 
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Here each U(l)i {i — 1, 2) corresponds to the phase rota- 
tion of *! or while U(l) mass and t/(l) sp in correspond 
to the overall and relative phase rotations, defined by 

C/(l)mass : *i -)■ *ie M , *2 -> * 2 e M , 
U(l) spin : * 1 ^* 1 e 4 ' 3 , tf 2 -> ^e"^, (12) 

whose currents are mass and pseudo-spin currents, re- 
spectively. Both the condensates \&i,\I r 2 are unchanged 
under the Z 2 action (a = f3 = 7r) inside J7(l) mass x 
t^(l)spin m Eq. (fT2"j). and therefore this Z 2 has to be re- 
moved, as the denominator of Eq. (111]) . 

In what follows, we call the phase cycles for \&i and 
*2 the (1,0)- and (0, l)-cycles, respectively. 



B. Vortex configuration 

Since the first homotopy group of order parameter 
space is 



7Ti(T 2 ) = Z©Z, 



(13) 
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it allows two kinds of winding numbers. We refer a vor- 
tex winding around (1, 0)[(0, l)]-cycle once as a (1,0)- 
vortex [(0, l)-vortex], which is the most fundamental vor- 
tex. When one travels around a (1,0) [(0, l)]-vortex, the 
phase of ^i(^> 2 ) rotates by 2ir with the phase of the 
other component kept constant. On the other hand, in 
terms of U(l) mass and U(l) spin in Eq. (TT2]), U(l) mass is 
rotated by ir and U (l) sp i n is rotated by +ir (— if) with cir- 
cling around a (1, 0)[(0, l)]-vortex. Since they have a half 
winding of ?7(l) mass , they are often called half- quantized 
vortices. 

Vortices winding around both components by 2tt are 
denoted by (1,1) and have unit winding in U(l) mass . 



They are called integer vortices, if the core is not sep- 
arated into (1,0) and (0, 1) vortices. More generally we 
refer a configuration which winds (l,0)-cycle m times 
and (0, l)-cycle n times as an (m, n)-vortex, whose wave 
function is denoted as vjo" 1 '") for i-th component. 

The vortex configuration can be obtained by solving 
Eqs. ([6]) and ((7|). Let us make an ansatz for an axially 
symmetric (1, 0)-vortex configuration 



(i,o) _ 
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where r and 9 are the polar coordinates. 



The profile functions /(i,o) and ft(i,o) are determined 
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with the prime denoting a differentiation with respect to 
r. We solve these equations with the boundary conditions 

(/(i,o)jV,0)) -> (M) as r->oo, (17) 
(/(i,o),ft'(i,o)) "> (0,0) as r^0. (18) 

From these equations, asymptotic behaviors of the pro- 
file functions /(i,o) and ft(i,o) at large distance can be 
obtained as 



/(i,o)(r) = 1- 
h(i,o)(r) = 1 + 
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where we have introduced the effective mass parameters 

+ 4(^152-^2512) _ 4(/x 2 5i - M1312) , 01 x 
^1 = TT2 . % = — is ■ ( 21 ) 
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The stability condition Eq. ([9]) of the ground state en- 
sures that rji > 0, while r]^ changes its sign with g i2 . 
Similarly, we make an ansatz for the (0, l)-vortex 



- «iV,i)W: = «ae w / (0l i)(r). (22) 

The equations for /(o,i), ft(o,i) can be obtained by just 
replacing the indices as 1 f-> 2 and (1,0) f-> (0,1) in 
Eqs. (fl"5|) and (fT6|) . Then the asymptotic behaviors are 
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Again, 77^ i s always positive while sign of r/ 2 depends on 

3l2- 

As vortices in a scalar BEC, the tension (energy per 
unit length) of (1,0) and (0,1) vortices logarithmically 
diverges as 



wh 2 v 2 L 
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respectively, with L and £ being a long and short distance 
cutoff, respectively. This divergent behavior comes from 
the kinetic term in the GP energy functional Eq. ([3]). 

Some numerical solutions of the single vortex configu- 
rations are shown in Fig. [TJ A universal feature of con- 
figuration is that ft(i ! o) (the profile function of unwinding 
field) at the vortex center is concave for 312 < and con- 
vex for the (?i2 > [3l| . This can be understood from 
the atom-atom interaction g\ 2 ; in the presence of the vor- 
tex profile for "fi as a background, ^> 2 feels the potential 
312 |^i 1 2 and it tends to be trapped in the vortex center 
for the repulsive interaction g\ 2 > and to be exclu- 
sive from the vortex center for the attractive interaction 
312 < 0. 



C. Inter- vortex forces 

It is expected that the interactions between (1,0)- and 
(0, l)-vortices are determined by the coupling (712-term. 



4 



-20 







--^.4- 





\ \ - 

U- 


1 
1 1 


l]p- 




o.y 




0.6- 




o.4 ; 




0.2 


gu = -0.3 



-10 



10 



20 



(a) 



1.4 




1.2 








Vs 

06 




o.k 




oi 







(b) 



1.4 






1.2 






1.1) 






/ 






_ -D.8 


\ 




\o.6 






OA 






0.2 




gi2 = 0.3 



-20 -10 10 20 

(c) 

FIG. 1: Single vortex configurations Q^il 2 (solid line) 
and [ \I> 2 1 2 (broken line)) on a cross section. The field 
winds once so that |>i?i| goes to zero at the vortex cen- 
ter while ^2 does not touch zero anywhere but it can have 
nonzero amplitude at the vortex center. The parameters are 
(h, gi,g 2 , pii,M2, mi, m 2 ) = (1,1,1,1,1,1,1) and (a) gi 2 = 
-0.3, (b) 5 i2 = and (c) g 12 = 0.3. 



When gi2 is zero, they are decoupled in Eqs. ((4]) and 
(O, so (1, 0) and (0, 1) vortices do not interact. Here, we 
calculate the asymptotic interactions between well sep- 
arated (1,0)- and (0, 1)- vortices using their asymptotic 
profile functions obtained in the last subsection. Let us 
place the (1, 0)- and (0, l)-vortices at (a:, y) — (R, 0) and 
{x 7 y) = (— i?,0), respectively as in Fig. O We use the 
polar coordinates (r, 9) with the origin (x, y) = (0,0). 
We further express (rnm, #(i,o)) an d ( r (o,i), ^(o,i)) as the 
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FIG. 2: Configuration of (l,0)-vortex and (0, l)-vortex. 



polar coordinates with the origins at the (1, 0) and (0, 1) 
vortex centers (R, 0) and (— R, 0), respectively. Then the 
following relations hold among three polar coordinates 
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with i = (1,0), (0,1), the minus sign for i = (1,0) and 
the plus sign for i — (0, 1). With these coordinates, the 
(1, 0)- and (0, l)-vortex configurations (^J 1 , $2 ) an d 
(\E , i°' 1 \ \E , 2°' 1 ^) can be expressed as 
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Let us now calculate the interaction between (1, 0)-vortex and (0, l)-vortex. We first make the standard Abrikosov 
ansatz 



*< M) (r,0) = v^VY'"'*?' 1 ' si wi [ 1 



for the total configuration. 
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Then the interaction potential is obtained by subtracting 
two individual vortex energies from the total energy as 



(32) 
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where it has two contributions: one from the kinetic 
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By using the asymptotic properties given in Eqs. (119[) . 
201), O, (El, (0Q1) and (JSJ, we find 
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where we have used (Vfl^.o)) 2 = r (i 2 o)> (V^(o.i)) 2 = 
r (oi) an< ^ nave taken terms up to 0(r~ 4 ). It is impor- 
tant to see that the leading terms of the order (D(r~ 2 ) 
have been canceled out in the subtraction. Therefore, 
the dominant contribution to the interaction potential is 
of the order 0(r~ 4 ). Similarly, we find the terms of the 
order 0(r~ 4 ) in the potential energy 
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Plugging these into Eq. (|3"2")) , we get 
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where £ stands for a short distance cut-off and we have 
used R >• £. The detailed calculation of Eq. (|35l is de- 
scribed in Appendix. Here, the terms independent of R 
have been ignored. The factor 1 /R 2 is a striking feature 
which is absent in the scalar BEC or scalar superfluids. 
Note that the chemical potentials \x\ and fi 2 do not ap- 
pear in the final result (|35l) . The asymptotic force be- 
tween the two vortices is obtained by differentiating the 
potential by their distance 2R, as Fni\(R) — -- '" "•" 



2dR 



F(IA){R) 



1 



4miTO 2 (gi32 - 3i 2 ) ^ 3 V i 



log 



R 



.(36) 



We have found that the interaction is attractive for g\ 2 < 
0, repulsive for 1712 > and vanishes for g± 2 = 0. 

Note that the asymptotic interaction is independent of 
the sign of the vortex winding number e ±l6 , because the 



interaction between the two condensates is mediated only 
through their amplitudes as (712 1 ^1 1 2 1 ^2 1 2 - In fact, the 
interaction potential ^ 1 , — 1 ) between (1,0) and (0,-1) 
vortices are exactly the same as Un iy It is easy to verify 
that the following relation holds 



^(1,1) = ^(1,-1) = ^(-1,1) = ^(-1,-1)- 



(37) 



This is because 0(i,o) and @(o,i) are decoupled in the 
Abrikosov ansatz in Eqs. (f3T)]) and (f3"Tj) . 

The potential (f35)) should be compared with the po- 
tential C/(i±i,o) between (1,0) and (±1,0) vortices. To 
see it, we make the ordinary Abrikosov ansatz 
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Note here that we have taken terms of the order unity. A 
leading order contribution to the interaction comes from 
the kinetic term of Vfi which is of order 0(r~ 2 ). On the 
other hand, the kinetic energy of Vf^ and the potential 
energy contributions start from the order 0(r~ 4 ), so we 
omit them. The interaction potential is then given by 
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where L is an infrared cut-off parameter; see Appendix 
for the details. Unlike the case of the leading term in 
the potential ((35)) between (1,0) and (0,1) vortices, the 
potential Un±i t o) depends on the chemical potential. We 
also note that it depends on the infrared cutoff L but not 
on the ultraviolet cutoff £. 



The inter-vortex force -F(i±i,o) 
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where — > denotes the large volume limit L — > 00. This 
l/R force is well known for vortices in the scalar BEC, 
scalar superfluids and the XY model, and global vortices 
in relativistic field theories [l|| [l8[ . The correction term 
for finite volume L can be found in the second term in the 
brace in the first line. This term might not be so familiar 
but has been obtained previously in [22j for global non- 
Abelian vortices in QCD. 

In the same way, the interaction potential between 
(0, 1)- and (0, ±l)-vortices are given by 



U {0A±1) (R)=± 



mgi2)h 2 -K R 2 

log 
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III. NUMERICAL ANALYSIS 



Let us numerically verify the interaction potential an- 
alytically obtained in Eqs. (|35|) . For simplicity, we con- 
sider a special case mi = 777,2 = ffi, gi = g 2 = g and 
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FIG. 3: The inter-vortex potential U {1>1) (R) for m_ = (0.5,0.6,0.7,0.8,0.9, 1, 1.2, 1.4, 1.6, 1.8,2) with m H 
asymptotic inter- vortex forces (Abrikosov ansatz) which are analytically obtained in Eq. (|35[). 
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= 1. Solid lines are 



/ii = fj& = fi. Then the asymptotic behaviors of the 
profile functions in Eqs. (fl"9j) and ([20)) are rewritten as 
follows, 



fi + hi = 2- 



(42) 
(43) 



for i = (1, 0) and (0, 1), where the mass parameters 
and m_ are defined by 



4m^i 



512 
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The inverse numbers of m + and m_ give the healing 
lengths associated with the mass component fi + hi and 
the spin component fi — hi, respectively. The inter- vortex 
potential Eq. (|35|) is then expressed as 



= „ 2 „2 7^ lQ g- 



(45) 



where = H 2 /j,/m(g + g 12 ) is defined in Eq. (|10l) . 

To obtain the inter-vortex potential numerically, we 
use a sort of the imaginary time propagation of the GP 



equation as 
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V 2 - /i + 9 1 * 2 (x) 1 2 + ,g 12 1 * ! (x) 1 2 ^ * 2 (x, r) 

= - J D 2 (x)a T * 2 (x,r),(47) 



where t is the imaginary time and D\ and Di are pos- 
itive coefficients. While Di is set to be a constant in 
the usual imaginary time propagation, we consider the 
coefficient Di = Di(x) with space-coordinate depen- 
dence. An advantage of using -Di(x) is that one can 
effectively fix the position of vortex (the zeros of 
during the numerical calculation, if one choose Di (x) ap- 
propriately fn order to attain this, we choose a function 
Dj(x) = AV 2 log (|x — ai\ 2 + e 2 ) + c where a, stands for 
the i-th vortex position and A and c are positive con- 
stants. The value of A is taken as an extremely large 
value to fix the profile of the wave function only near the 
vortex cores. Also, e should be sufficiently small. We 
chose A = 80000, e = 0.01 and c = 0.1 in our numeri- 
cal computation. We take the Abrikosov ansatz given in 
Eqs. (f3"0"|) and (|3"Tj) as the initial condition at r = and 
minimize the energy under the imaginary time evolution. 
After the solutions converge sufficiently, we calculate the 
interaction energy Eq. (f3"2"l) . 

Throughout our numerical computation below, we will 
set m/h 2 = 1 and v 2 = 1. Then we regard m + and m_ as 
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independent parameters of the GP equations and perform 
the numerical calculation by varying them. Remember 
that m+ < to_ corresponds to g±2 < (attractive force) 
whereas m + > m_ corresponds to gi2 > (repulsive 
force). No net interaction exists accidentally when m + = 
m_. 

The result is shown in Fig. [31 We compare the inter- 
vortex potential obtained numerically and the one ob- 
tained analytically. As can be seen, the analytic results 
reproduce the numerical results quite well. We have only 
one fitting parameter 2£ which is the the short-range cut- 
off. The values 2£ for various choice of m_ for fixed 
m + = 1 arc shown in Fig. |4j We find the linear depen- 
dence of the short-range cut-off £ on the healing length 
l/m_. 

2f 
0.9^ 

0.8 ; 
0.7 ; 
0.6; 

0.5 \ 
OA 



0.8 1.0 1.2 1.4 1.6 1. 



2.0 



FIG. 4: The dependence of the healing length 2£ on l/m_ 
with fixed m + — 1. 



IV. SUMMARY AND DISCUSSION 

We have studied the asymptotic interaction between 
half-quantized vortices, i.e., (1,0) and (0,±1) vortices 
winding around "Ji and ^2, respectively in the two- 
component BEC. Since the two components interact only 
through the density, the (l,0)-vortex does not directly 
experience the circulation of the (0, ±l)-vortex so that 
the result does not depend on the signature of the wind- 
ing number. The leading order of the force between them 
is found to be ~ [log(i?/£) — l/2]/R 3 in contrast to the 
one between the same kind of vortices ~ 1 / R, which is 
also well known as the force between vortices in scalar 
BEC, scalar superfluid and the XY model, and global 
vortices in relativistic field theories. We have first de- 
rived it analytically using the Abrikosov ansatz and the 
asymptotic profile functions of (1,0)- and (0, l)-vortices. 
We have then confirmed it numerically with using the 
extended imaginary time method for the GP equations. 
We have found that the short-range cut-off parameter £ 
of the vortex interaction linearly depends on the healing 
length l/m_. 

Our results suggest a bound state of (1, 0) and (0, ±1) 
vortices for c/12 < 0. While a set of (1,0) and (0,1) is 



expected to form a stable integer (mass) vortex (1,1), 
it is a nontrivial question if (1,0) and (0,-1) vortices 
form a bound state, which should be called a (pseudo-) 
spin vortex. Also, one expects no stable bound states for 
.912 > 0. Although there should be instabilities for large 
separation at least, it does not exclude a possibility of 
a metastable bound state at short distance. To address 
these questions, we need to know a short range inter- 
action or stability analysis of the bound states, which 
remains as a future problem. 

Multiple vortices will constitute a vortex lattice in ex- 
periments of multicomponent BECs under the rotation 
32]. Ample phase diagram of the vortex lattices was 
predicted in [33| and was numerically obtained in 0, l3~ij . 
Since we have obtained the analytic expression of the 
inter-vortex forces, we expect to explain the vortex phase 
diagram analytically. Especially, we expect that the dif- 
ference of the force [log(i?/£) - l/2]/i? 3 between (1,0) 
and (0,1) vortices and the one 1/R between the same 
kind of vortices will determine it. 

Our method should be extended to spinor BECs, which 
remains as an interesting future problem. On the other 
hand multicomponent systems in relativistic field theo- 
ries are common in QCD such as the linear sigma model 
for the chiral phase transition and the Landau-Ginzburg 
model for color superconductors at high baryon density 
23]. In these models, order parameters are matrices as 
in superfluid 3 He rather than vectors, and consequently 
there exist non-Abelian vortices [35| : non-Abelian global 
vortices in the chiral phase transition 211] and non- 
Abelian semi-superfluid vortices in color superconductors 
at high baryon density [23[ . Inter- vortex forces have been 
calculated at leading order for non-Abelian global vor- 
tices i|22[ and non-Abelian semi-superfluid vortices [24| . 
see [25| for a review. Calculation in the present paper 
will give the next leading order [log(i?/£) — l/2]/i? 3 to 
them. Especially the force cosa/i? between non-Abelian 
global vortices at the leading order vanishes for a partic- 
ular choice (a = ±7r) of internal orientations of vortices 
[22| . and therefore the next leading order term propor- 
tional to [log(i?/£) — 1/2]/ R 3 becomes a dominant contri- 
bution. An extension of our results to these cases should 
be important to consider a possibility of vortex lattice 
phases in heavy-ion collisions or in a neutron star core, 
as in two-component BEC @, HH, [34| . 
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FIG. 5: The integral region to calculate Eq. (|A2 



Appendix A: Derivation of Eqs. (|35[) and (|39|> 



In Eq. (|35|) , we have to evaluate the integral 



/ = I d 2 x 



An 



r 2 r 2 / /!„ ' 

r (i,o) r (o,i) ^ A ° 

.2 D 2 . 



i? 4 - 2r^i? 2 cos 26». 



To this end we will use a formula 



dO- 



2tt 



A + Bcos 29 ^/A 2 - B 2 ' 



(Al) 
(A2) 

(A3) 



for A > To evaluate Eq. (|A2I) . we divide the integral 
region as shown in Fig. [SJ In addition to I\ and Ii-, 
we take into account the contributions I3 from the strip 
of width 2£. Since the integrand diverges at (x, y) = 
(±i?, 0), we introduce an ultraviolet cut-off £. Then, we 
will remove the small regions that includes the points of 
vortex positions. Hence, the total integral is written as 



cut off 



h 



2/3 



(A4) 



For R ^> £ the integral l\ and I2 is calculated as 
" ~Z 2irrdr 



h = 



h = 



(A5) 



27rrdr tt / R m ( $ \ . , 



where we have used Eq. (|A3[) . The remaining integral 



r 

d9— 

An 



(A7) 



dr 

with <C 1 and <5 <gc 1 is evaluated as follows. Note 

that cos 29 < cos 26 = 1 — 26 2 H < 1 - J 2 and A > 

(r 2 - i? 2 ) 2 + 2r 2 i? 2 <S 2 > 2r 2 R 2 6 2 . Thus, we have the 
following inequality 



„ T IT -25, i? + £ 7T £ 

< / 3 < log 



2R 2 6 2 



R-Z R 2 6 2 R 



(A8) 



Thus, for any 5, one can choose sufficiently small £, so 
that J3 becomes negligibly small. 
In summary, we get 

£(*f + o(I)). (A9) 

Next, we calculate the integration in Eq. ([3T)|) . 

r 2 -i? 2 



J = I d 2 x V6>(i, ) • V^(o.i) = I d 2 x — 



drd9- 



r{r 2 - R 2 ) 



R 4 - 2r 2 R 2 cos 26*' 



(1,0)' (0,1) 

(A10) 



By using Eq. (|A3I) . one can first perform the integration 
in 9 and then integrate with respect r as 



J = 



dr 



2nr{r 2 - R 2 ) 



= lim 

L— ¥00 



TF^W 

fl 2?rrdr 
KR 2 







L 2ixrdr 
> r 2 + i? 2 



, i 2 + i? 2 

lim 7r log = — . 

L^oo B 4i? 2 



(All) 
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